GEOMETRIC CONTROL METHODS FOR QUANTUM 
COMPUTATIONS 



ZAKARIA GIUNASHVILI 

Abstract. The applications of geometric control theory methods on Lie groups 
and homogeneous spaces to the theory of quantum computations are investi- 
gated. These methods are shown to be very useful for the problem of con- 
structing an universal set of gates for quantum computations: the well-known 
result that the set of all one-bit gates together with almost any one two-bit 
gate is universal is considered from the control theory viewpoint. 



1. Introduction 

There are three main directions in the theory of quantum computations: 

• Algebraic: the problem of constructing a universal set of gates 

• Physical (practical): the problem of constructing a physical device imple- 
menting a universal set of gates 

• Algorithmic: inventing algorithms for quantum devices to solve various 
computational problems effectively, provided that a universal set of gates 
is already constructed (at least theoretically). 

The first of these directions is pure algebraic (perhaps linear algebraic), since a 
universal set of elementary gates is such subset of the group of unitary transfor- 
mations of some finite-dimensional Hilbert space which generates the entire group. 
In the case of quantum computer there is a correspondence between computational 
tasks and unitary transformations, and the quantum computer program solving 
the problem is a composition of some elementary gates which equals to the unitary 
transformation corresponding to the given problem. Therefore, any conventional 
(as classical so quantum) computational process is discrete while the evolution of 
any physical system (as classical so quantum) is continuous and is described with 
the corresponding Hamiltonian or Schrodinger equation. Hence, it can be stated 
that during a conventional computational process, when we want to represent a 
given unitary transformation as a composition of some elementary gates, we have 
a mechanical system with Hamiltonian controlled by some external means (inputs) 
and we are able to implement the elementary gates by these inputs. The compu- 
tational process is considered as discrete because we are interested not with the 
continuous evolution of the system but with the transformations at some finite mo- 
ments of time corresponding to the elementary gates. Nowadays there are several 
works where the "discrete" approach to the computational process is replaced by 
"continuous", when the aim is to construct a given transformation of a system not 
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by means of discrete inputs generating the preset gates, but by continuous steering 
with some external parameter-dependent Hamiltonian (very interesting approach 
is developed in P^, [2], [3]). The geometric control theory is a tool which unifies 
these two approaches and the aim of this work is to demonstrate the methods that 
can be used to build a computational process (program) by parameter-dependent 
Hamiltonian and to investigate what can be achieved by these methods. 

2. Dynamical Polysystems and Control Group 

Definition 1 (Control Group, see 0]). Let I be any nonempty set. Consider the 
set of all finite sequences of type {{t\,i\) ■ ■ ■ (t m ,i m )) where (tk,ik) € R. X I, k = 
1, . . . , to. Introduce the following reduction rules: 

(1) any term of the type (0, £) is removed from the sequence 

(2) if ik = ifc+i = i then the pair of consecutive terms (tk,ik)(tk+i,ik+i) is 
replaced by the single term (tj. + tk+i,i) 

It is clear that after a finite number of steps from any sequence can be obtained an 
irreducible sequence (which can be empty and we denote such one by 0). We denote 
the set of all irreducible sequences by C(I) and call the control group over the set I . 

Remark 1. Another point of view is to consider the control group as a quotient 
set of the set of all sequences {((ii, £].)••• (t m , i m ))} under the following equivalence 
relation: two sequences are equivalent if after reduction they give one and the same 
irreducible sequence. 

The following operations on the set C(I) justify the usage of term group. For 
any two sequences Si and s 2 from C(I) let (S1S2) be the sequence obtained by their 
concatenation and let the product si ■ S2 be the irreducible sequence obtained by 
the reduction of (sis?). The operation (si, S2) 1— > s% ■ S2 together with the neutral 
element (the empty sequence) defines a group structure on C(I). 

Remark 2. The group C(I) is commutative only when I is a one-element set. 

The elements of C(I) are called controls. For any control s = ((ti, oq) • ■ ■ (i„, a n )) 
its inverse element in the group C(I) is the control s _1 = ((— t n , a n ) ■ ■ ■ (— ti, ai)). A 
control ((ti, ai) • • • (t n , a n )) is said to be positive [negative) if ti > (U < 0), Vie 
{1, . . . ,n}. Let us denote the set of all positive (negative) controls by C + {I) (C~(I)). 
In general C(I) ^ C + (I) U C~(I) because there are many "mixed" type controls in 
C(I). It is clear that C + (I) and C~(I) are semigroups. 

For any A G R and s = ((ti,ai) ■ ■ ■ (t n ,a n )) 6 C(I) define a control A • s as 
A • s — ((Aii, o>i) • • • (Ai„, a n )). Any subset {s%, . . . , s p } C C(I) defines a mapping 

$ sx ... Sp : W — ► C(I), $ ai ... Sp (Ax,...,Ap) = A lSl • \ 2 s 2 ---\ps p 

After this we can introduce a topology on C{I) as the strongest topology for which 
all the maps of type ^^...s are continuous. Analogically can be introduced a 
differential structure on C(I): a mapping f : C(I) — ► M is said to be smooth if the 
mapping / o $ Sl ... s : R p — > M is smooth for any $ Sl ... s . The group operations 
in C{I) are continuous and smooth with respect to the above defined topology and 
smooth structure. 

Definition 2. For a smooth manifold M a dynamical polysystem on M controlled 
by C(I) is a smooth action of the group C(I) on M : 

C(J) x M 3 (s,m) i-> sme M 
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Definition 3. The set C(I)x — {sx \ s G C(I)} is called the orbit of the point 
x G M. The set C + (I)x = {sx | s € C + (I)} is called the positive orbit of x and 
C~(I)x = {sx | s G C~(J)} - the negative orbit of x. 

In this case any fixed element a G / gives a one-parameter subgroup C a = 
{(t, a) t G K} C C(I) which induces a one-parameter group of diffeomorphisms of 
M: 

ift : M — > M, t G K, tp%(x) = (t, a)x, x G M 

The latter itself corresponds to some smooth vector field X a on M. The family of 
vector fields {X a \ a G /} is called the infinitesimal transformations of the dynam- 
ical polysystem. Conversely, if we have a family of vector fields F on M indexed 
by the elements of some set /, then we obviously have a dynamical polysystem on 
M controlled by C(I). In such cases we denote the control group together with 
the corresponding dynamical polysystem by C(F). Thus, we can state that there 
is a one-to-one correspondence between the dynamical polysystems on M and the 
families of smooth vector fields on M. 

3. The Structure of an Orbit of Dynamical Polysystem 

Let F be a family of vector fields on the manifold M . 

Definition 4. The closure of the family of vector fields F is the family of vector 
fields defined as 

A(F) = {sX | s G C(F), X G F} 
In other words it is the minimal set of vector fields containing F and invariant 
under the action of the control group C(F). 

The dimension of the dynamical polysystem C(F) at a point x G M is the 
dimension of the subspace of the tangent space T X (M) generated by the family of 
vector fields A(i ? ). The main result about the structure of an orbit of dynamical 
polysystem is based upon the following 

Lemma 1 (see 5 ). For any point x G M the dimension of the dynamical polysys- 
temC(F) is constant through the orbit C{F)x. If m is the dimension of the dynami- 
cal polysystem C(F) at the point x then there exists a set of controls {s%, . . . , s m } G 
C{F) and a set of vector fields {Xi, . . . , X m } C F such that the mapping 

F ■ E m — ► C(F)X, ...,t m ) = (Smells- 1 ■ • -Sl^sr 1 ) 

where {ip\ \ t G M} , i G {1, . . . , to} is the one-parameter group of diffeomorphisms 
corresponding to Xi, is a local diffeomorphism at the point (0, ... ,0) G R m . 

Let exp(F) denote the group of diffeomorphisms of M generated by the flows of 
the elements of -F: 

exp(F) = {exp(tiXi) • • ■ e^p(t n X n ) | t 4 G M, X. L G F, n G N} 

and C(F) be the minimal Lie subalgebra of T{M) containing F and invariant under 
the action of the group exp(F) . 

From the above lemma follows the following 

Theorem 2 (see [S]). For any point x in the orbit C(F)xq, xq G M , the mapping 
<f) x : M m — > C(F)xq gives a local coordinate system on C(F)xq. The topology and 
the differential structure induced on the orbit by these coordinate system coincide 
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with the topology and the differential structure induced by the action of the group 
C(F). 

and 

Theorem 3 (Orbit Theorem, Nagano-Sussmann). For any point xq G M the orbit 
C(F)xq = XQ is a connected immersed submanifold of M and T X0 (O X0 ) — C(F) Xo . 

4. Control System and Accessibility Problem 

For a smooth manifold M let T(M) be the space of smooth vector fields on M. 
We shall generalize the notion of control group in the following way. Let U be some 
nonempty set (usually assumed to be a finite-dimensional smooth manifold). We 
denote by V(U) the set of all paths in U. At this point we do not specify what 
type of paths we consider 

V{U) = {u : [0, T] — > U | T > 0} 

The set V{U) is a semigroup under the following standard operation: for u\ : 
[0, Ti] — > U and u 2 : [0, T 2 ] — ► U let u 2 u 1 : [0, T x + T 2 ] — ► U be 



(«2«l)(*) = 



ui(t), te[0,Ti]; 
u 2 {t), t G [Ti, Ti + T 2 ]. 

A control system on the manifold M is the following data 

• a smooth map A : C7 — > JT(M); 

• a subspace C(f7) of 7- > (C/) (at this stage we do not specify the structure of 
C(U), but it is assumed to be a subsemigroup of V{U)). 

The elements of C(U) are called controls and any u G C(U), u : [0, T] — > t/ defines 
a dynamical system 

(1) m = X u (t) m(t) , m(t)eAf, te[0,T] 

where X" denotes the time-dependent vector field A"(£) = (Xozt)(t) and X"(t) m ( t ) 
is the value of this vector field at the point m{t) € M. 

Remark 3. The solution o/Q is defined as the absolutely continuous map m : 
[0, T] — ► M such that the equation is satisfied for almost every t G [0, T]. 

For any given control u : [0, T] — ► U and an initial value m(0) = mo, we denote 
by j(mo, u, t), t G [0, T] the solution (if exists) of^ If for some t G [0, T] and x G M 
we have that 7(7710, h, t) — x, then we state that the control u steers the state mo 
into the state x in time t, and the state x is reachable from mo in time i. We denote 
by 1Z(xq, t) the set of points of M reachable from the point xo G M in time i El. 
Traditionally also the following two sets are the subjects of investigation 

it(xo,t) = [J 1Z(xq,t) - the set of points reachable from xq in time < t 

Q<T<t 

and 

1Z(xq) = TZ(xq, t) - the set of points reachable from xq 

0<T<OO 

The aim of geometric control theory is the investigation of the structures of the 
sets 7Z(xo,t), lZ{xQ,t) and TZ(xq). These problems are considered for various types 
of the space of control parameters U as well as the map X : U — > J-{M) and 
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the class of controls C(U). We assume that U = R m for some integer to > and 
consider the following three types of controls 

• C u - unrestricted controls, which consists of locally bounded and measurable 
mappings from [0, T] to R m , T > 0; 

• C r - restricted controls: the subset of C u consisting of the mappings with 
values in the cube {(xi, . . . , x m ) G R m |||xj| < 1, i = 1, . . . , to}; 

• Cb - "bang-bang" controls, consisting of the all piecewise constant maps 
from [0,T] to R m . 

It is clear that in the case of bang-bang controls the reachable set is independent of 
the structure of the map X : U — ► T(M) and depends only on the image of this 
map. The set reachable from a point to G M in time t > for piecewise constant 
controls is 

TZ(m, t) = {exp(iiXi) • • • exp(t k X k )m \ 

k 

E *« = *> U> 0, X l G Image(X), i=l,...,k, keN} 
i=i 

where exp(i£) denotes an element of the one-parameter group of diffeomorphisms 
associated with a vector field £. In the case when the set Image(X) is symmetric in 
the sense that £ G Image(A) — £ G Image(X), we have that the set reachable 
from any point to G M coincides with the orbit of the dynamical polysystem on 
M defined by the family of vector fields Image(X). So we can apply the classical 
Orbit Theorem (see Theorem EJ for the structure of the set lZ(m), m G M. 

Let Cie(S) be the minimal Lie subalgebra of J-(M) containing the subset S C 
T(M). It is clear that for any point x G M the space Cie(S) x C T X (M) is a 
subspace of the tangent space to the orbit of x. In some cases it coincides with the 
tangent space of the orbit. Particularly, the following theorem is one of the most 
important corollaries of the Orbit Theorem 

Theorem 4 (Rashevsky-Chow). If M is a connected smooth manifold and Cie(S) x = 
T X [M) for every point x G M , then the orbit of the point x under the action of the 
dynamical polysystem defined by S coincides with the entire manifold M . 

A set S G T{M) which satisfies the condition Cie{S) x = T X (M),V x G M is said 
to be completely nonholonomic or bracket- generating. As it was mentioned if the 
set S is symmetric then the orbits of the corresponding dynamical polysystem and 
the reachable sets coincide. Hence, we can conclude that for symmetric (S = —S) 
nonholonomic family of vector fields the corresponding control system is completely 
controllable. 

5. Two Notions of Controllability in the Case of Classical 

Hamiltonian System 

Sometimes it is useful to distinguish the following two notions of controllability 
(reachability) as for classical so for quantum mechanical systems. 

Definition 5 (Transformation-Controllability). We call a system transformation- 
controllable if every transformation of the system can be obtained by applying of 
some finite set of appropriate controls in a finite interval of time. 

Definition 6 (State-Controllability). We call a system state-controllable if for 

every pair of states of the system - s± and S2, there exists a finite set of controls 
that steers the system from the state s± into the state S2 in a finite time. 
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In the case of classical system, if we assume that the set of all transformations 
is the entire group of diffcomorphisms of a smooth manifold M - Diff(M), and 
the control system is the dynamical polysystcm defined by some family of vector 
fields F, then the transformation-controllability means that exp(F) = Diff(M); 
and the state-controllability means that for any point m G M: O m = M. These 
two notions are not equivalent neither in classical nor in quantum cases. To demon- 
strate this let us consider the following situation. Let M be a symplectic manifold 
and the dynamical polysystem F be the family of all Hamiltonian vector fields 
on M: F — ham(M). Let { , } denote the Poisson bracket corresponding to 
the symplectic structure on M. According to the Darboux theorem for any point 
m G M there exists a local coordinate system {pi, q{ \ i = 1, . . . ,n} such that 
{PiiQj} = This implies that locally we have the set of Hamiltonian vector 
fields |^7) | i = 1, . . . ,n| which forms a basis for T m (M). So we have that 

for any point m G M the set of Hamiltonian vector fields forms a basis of T m (M). 
Therefore, by the theorem of Rashevsky-Chow the family of all Hamiltonian vec- 
tor fields is completely nonholonomic and therefore the corresponding system is 
state-controllable. But the family ham(M) is not transfomation-controllable be- 
cause the corresponding family of 1-parameter flows is contained in the group of 
symplectic diffeomorphisms of the manifold M, and it is clear that this group is 
just a part of the entire group of diffeomorphisms of M. The classical problem of 
state-controllability for the family of all Hamiltonian vector fields becomes more 
complicated if we extend the notion of state and consider not only pure states (i.e. 
points), but also mixed states. Let us recall the formal definition of this notion 
which is rather well-known from statistical physics. The realization of the alge- 
bra of observables as the algebra C°°(M) dictates the realization of the space of 
states as the space of probability distributions, i.e., the linear positive functional 
p : C°°(M) — > M such that <j) > p(<j>) > and / p{p, q)dpdq = 1. Let us dc- 

M 

note the space of such functionals by Stat(M). In other words, a state p G Stat(M) 
is a positive distribution (generalized function) such that p(l) = 1. The pure state 
corresponding to a point x G M is the Dirac functional 5 X : 5 x {<f) = 4>(x). Some- 
times we shall use the scalar product notation (p, <fi) for p{4>). Any smooth map- 
ping of pure states / : M — ► M can be extended to / : Stat(M) — ► Stat(M) as 
(f(p), 4>) = (p, /*((/>)) where /* : C°°(M) — ► C°°(M) denotes the mapping dual 
to/. 

Remark 4. It is clear that the subspace Stat(M) in the space of all general- 
ized functions is invariant under such mappings, i.e., every positive distribution 
is mapped to a positive one and the property p(l) = 1 is preserved. 

After this we extend the problems of state-controllability and reachability to the 
space Stat(M). First of all let us notice that the pure-state-controllabilty does not 
imply the state-controllability in the extended state space. The simplest example 
is the following: let M be a simplectic manifold and u> be the symplectic form on it, 
let a — a ■ u> n , 2n = dim(M) be the volume form normalized so that J a — 1; take 

M 

the state £ which is a generalized function defined as (£, 4>) — J <j>-a, <fi £ C°°(M). 

M 

As it was mentioned, for the dynamical polysystem generated by the family of 
all Hamiltonian vector fields - Ham(M), the space of pure states is controllable 



GEOMETRIC CONTROL METHODS FOR QUANTUM COMPUTATIONS 



7 



but it is not true for the space Stat(M) because the state £ is invariant under 
the action of Hamiltonian flows. The reasonable classification of the entire space 
Stat(M) from the point of view of reachability under the dynamical polysystem 
defined by Ham(M) seems very difficult. As an example of application of Orbit 
Theorem we consider a simplified version of this problem. Particularly, we consider 
the reachability problem for the space of discrete mixed states. 
The space Stat(M) is a convex space: 

Pi,p 2 G Stat(M) c Pl + (1 - c)p 2 G Stat(M), V c G [0, 1] 

Therefore, for any finite set of points P C M and a function s : P — > M such that 
ip > and s ( x ) — 1) we can construct the distribution p{P,s) — s ( x )$x G 

Stat(M). Let us denote the subspace of all such states in Stat(M) by Stato(M). 
From the property f(S x ) = Sf( x ) follows that Stato(M) is an invariant subset of 
Stat(M). Now we consider the reachability problem for the dynamical polysystem 
defined by Ham(M) just for the space Stato(M). Particularly we describe the quo- 
tient space Stat (M)/ ~, where two states pi and p 2 from Stat (M) are equivalent 
iff there exists a finite set of Hamiltonians {Hi, . . . , H k } C C°°(M) and a set of 
time intervals {t\, . . . ,t k } C R such that exp(t 1 Vf/ 1 ) • • • exp(t k Vn k )pi = p 2 , where 
denotes the Hamiltonian vector field corresponding to the function Hi. 
For any p = p(P, s) G Stat(M) we call the set of points P C M the support of 
the distribution p and denote it by Supp(p). For any c G R we call the number of 
elements in the set s (c) the multiplicity of c and denote this number by m(c). 
For a given p(P,s) G Stat(M), the set of all pairs {(c, m(c)) | c G Image(s)} we 
call the spectrum of the distribution p and denote this set by Spec(p). It is easy 
to verify that for any diffeomorphism (or even bijection) / : M — > M we have 
that f(p(P, s)) — p(f(P), s o /). The latter implies that if two states p\ = p(Pi,si) 
and p 2 — p(P2,S2) are equivalent then Spec(pi) = Spec{p2)- It turns out that the 
equality of spectra is also a sufficient condition for the equivalency of two states 
in Stato(M). Before we start to prove the sufficiency let us proof the following 
auxiliary 

Lemma 5. Let Xi, . . . , X k be tangent vectors at k different points {mi, . . . , m^} C 
M. There exists such Hamiltonian H G C°°(M) that for the corresponding Hamil- 
tonian vector field Vh we have: Vr (mi) — X i} i = 1, . . . , k. 

Proof. For each point rrii take an open neighborhood Ui with canonical coordinate 
system {pj,q l j \ j = 1, . . . ,k} so that UiCiUi — when i ^ I. Since the system of 
all Hamiltonian vector fields is completely nonholonomic, for each Xi there exists a 
function f, G C oc (M) such that V ft (m,) = X t and /;(M \ Uf) = 0. It is clear that 
the function fi does not "affect" the function /; inside the neighborhood Ui when 

k 

i ^ I. Therefore the function H = Y^, fi is the desired Hamiltonian. □ 

i=l 

Theorem 6. Two states pi and P2 from Stato(M) are equivalent if and only if 
Spec(p 1 ) = Spec(p 2 ). 

Proof. Assume that the spectra of p\ — p(P\,si) and p 2 = p(P2, S2) are equal. For 
pi consider the ordered set (xi, . . . ,x n ) where n is the number of elements in Pi, 
Xi G Pi, i = l,...,n, and si(xi) < si(xi+i). In the same way we construct the 
ordered set (yi, . . . ,y n ) for p2- From Spec(pi) = Spec(p2) follows that si(xi) = 
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S2{yi), i = 1, ... ,71. Thus the problem is reduced to the existence of such element 
g of a Hamiltonian flow that g{xi) = g(yi), i — 1, . . . , n. 

For n G N let M n be the subset of M n — M x • • ■ x M defined as 

M" = {(pi, ...,p n ) \ Vi+ Pj when i ^ j, i,j = 1, . . . , n} 

M n is an open subset of M n and therefore it is a submanifold. The set M n \ M n 
is a union of submanifolds M™- = {(pi, . ■ ■ ,p n ) | Pi = Pj} , i ^ j- It is clear that 
dim(M™,) = (n — 1) • dim(M) and hence if dim(M) > 2 and M is connected then 
M" is also connected. But recall that M is assumed to be a symplectic manifold, 
so dim(M) > 2. Thus we have a connected manifold M n and the control group 
generated by the Hamiltonian vector fields on M acts on M n as 

g m . _> M n , g m {x u . . . , x n ) = (g( Xl ), . . . , g{x n )) 

where g : M — > M is an element of some Hamiltonian flow. The question is: is this 
action transitive or not? For a given Hamiltonian H G C°°(M) the action of the 
corresponding flow on M n is (afi, . . . ,x n ) i— » (exp(t_ff)a;i . . . ,exp(tH)x n ), ( 6 1: 
therefore, the corresponding vector field on M n is VJj = (Vh, ■ ■ ■ , Vh)- That is, 
for any point (x\, . . . , x n ) E M n the value of this vector field at this point is 
(Vh(xi), . . . , Vff(x n )). According to the Lemma[S]for any tangent vector (Xi, . . . , X n ) 
at the point (xi, . . . , x n ) there is a Hamiltonian H such that Vh(xi), ■ ■ ■ , Vh(x u ) — 
(Xi, . . . ,X n ). This implies that the dynamical polysystem on M n generated by 
the diagonal actions of the Hamiltonian flows on M is completely nonholonomic. 
Hence, any two states from Stato(M) with equal spectra are equivalent. □ 

Extending this to the entire set of states Stat(M) just by direct analogy, we 
can see that if two distributions u)\ and u>2 are equivalent then there exists such 
diffeomorphism g : M — > M that g(Supp(u>±)) = g(Supp(uj2)) ■ The problem of 
existence of such diffeomorphism that maps one given closed set to another given 
closed set is itself quite complicated. 

6. Control Systems on Lie Groups: Homogeneous and Affine Cases 

In the case of control systems on manifolds the main difficulty for controllabil- 
ity and accessibility problems comes from the fact that even for a small number 
of control parameters the control group is a subgroup of an infinite-dimensional 
Lie group. For control systems on finite-dimensional Lie groups when we restrict 
ourselves with right-invariant vector fields the problem becomes easier sometimes 
and can be solved completely. In this section we review the cases of homogeneous 
and affine controls which are discussed in the classical paper of Jurdjevich and 
Sussmann (see [H]). We start from the definitions of basic notions. 

Let G be a finite-dimensional Lie group and g be its Lie algebra. According to 
the general definition (see Section 0} we need a set of control parameters U which 
is assumed to be M. m for some positive integer m; a set of controls C(U) which is a 
subset of the set of mappings u : R + — > R m and usually is the set of continuous, 
piecewise-continuous, piecewise-constant or smooth maps; and a continuous map 
H : W n — ► T(G). In this case we consider such H that takes values in the subspace 
of right-invariant vector fields on G. Thus, we can define H as a continuous map 
H : W n — ► g and call it the Hamiltonian (with hope that it will not cause any 
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confusion). These data defines a family of control-dependent differential equations 
on G 

(2) g{t) = H(u(t))g(t), g(t) € G, g(0) = g , t > 

where u : R + — ► M m runs through the space of controls C(R m ). The right-hand 
side of the equation is a time-dependent right-invariant vector field on G defined 
by the Lie algebra elements H(u(t)), t G K + . As in general case the solution of the 
equation [21 is absolutely continuous map g : R + — ► G that satisfies this equation 
for almost every t > 0. For any given control u and initial state go we denote by 
7((?o> u, t), t G M. + the solution of the equation|2 

If for some 7(50, u,t) — g for some u G C(K m ) and £ > then we say that the 
control u steers the state go into the state g in time t and the state g is reachable 
from go in time t. As in general case, we denote by: 

• lZ(go,t) the set of all reachable in time t elements g G G; 

• 7l(go,T) = U lZ(go, t) - the set reachable in time < T; 

a<t<T 

• IZ(go) = U £) ~ the set reachable from <7 - 

0<t<oo 

From the fact that the system is right-invariant easily follows the following proper- 
ties of the reachable sets 

TZ(g,t)=n(l,t)-g, ti(g,T)=K(l,T)-g, K(g)=K(l)-g 

Therefore, for the right-invariant control systems on Lie groups it is sufficient to 
study the set reachable from 1. 

In the above-mentioned work of Jurdjevich and Sussmann ([S]) is considered the 
case when the Hamiltonian H : W n — > g is an affine function 

m 

(3) H(ui, u m ) = Xo + ^m ■ Xi, for fixed {X Q , X m } C L 

i=l 

We make a brief review of the results and methods for such control systems, because 
it turns out that they are useful not only for this specific class of systems but for 
other classes of control systems. 

Let us introduce the following three types of control parameters 

• C u - unrestricted controls: locally bounded and measurable mappings from 

R+ — > M m ; 

• C r - restricted controls: a subset of C u consisting of the mappings with 
values in the cube {(xi, . . . , x m ) G R™ | \xi\ < 1, i = 1, . . . , m}; 

• Cb - "bang-bang" controls: the subset of all piecewise constant mappings 
(in fact it suffices if \xi\ = 1, i = 1, . . . ,m); 

For these three classes of controls we use the notation C without subscript when 
some statement is formulated for each of them. 

The following theorem from the classical theory of Lie groups is one of the 
important tool for studying control systems on Lie groups. 

Theorem 7 (Yamabe's Theorem, see [J]). Let G be a Lie group and let H be a 
path- connected subgroup of G. Then H is a Lie subgroup ofG. 

The following theorem states that for affine Hamiltonian (seeOJ) the solution of 
the controls system exists and moreover is complete. 
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Theorem 8 (see For every control u G C and go € G there exists a unique 
solution g(t), t G K + o/0 defined for every t G [0, oo] suc/i that g(0) = go. 

Let 2 be the Lie subalgebra of L generated by {Xo, . . . ,X m } and G be the 
corresponding connected Lie subgroup of G. 

Theorem 9 (see ©). The set 72(1) is a semi-group in G and if it is a group then 
it coincides with G. 

Proof. Let g\ — 7(l,ui,ii) and 172 = 7(1,^2,^2) be two elements of G reachable 
from 1 applying the controls u\ and 112- Define a new control v as 



v(t) 



Ui(t) when r G [0, t{\ 

U2(t — ij.) when r > t\ 

Since the system is right-invariant we have that 7(1, v, t\ + t^) = giQ\- 

Now assume that 72(1) is a group. As the space of controls is path-connected, 
the subset 72.(1) is also path-connected and therefore, according to the Yamabe's 
theorem (see [7J) it is a Lie subgroup of G. Since 72.(1) G G, we have that its Lie 
algebra V is a subalgebra of g. Conversely, consider any element v(a±, . . . , a m ) — 

m 

Xq + aiXi 6 j, a, = ±1. By definition of 72(1) we have that exp(iu) G 72(1) 

i=l 

for t > 0, but since 72(1) is a group we have that exp(to) G 72(1) for t < too. 
The elements v(a±, . . . , a m ), where (a\, . . . , a m ) G K m , are generators of g which 
implies that g C V. Since the groups 72(1) and G are connected we have that 
72(1) = G. □ 

An affine control system X = {Xo,X\, . . . , X m ) is said to be homogeneous if 
Xo = 0. In other words the Hamiltonian depends on control parameters linearly: 

m 

H(m, . . . ,u m ) = J2 u i x i- 

i=l 

Theorem 10 (see jH]). For homogeneous right-invariant control system on a Lie 
group G the reachable set 72(1), for any of the three classes of controls, is the 
connected Lie subgroup of G corresponding to the Lie algebra generated by X. If 
we use the class of unrestricted controls C u then for each T > we have that 
72(1, T) = 72(1). 

Proof. Since the set 72(1) is a semigroup, according to the previous theorem it is 
sufficient to show that if g G 72(1) then g _1 G 72(1). Assume that g — 7(1, u,t). 
Consider a control v defined as 



v(t) 




t) when r G [0, t] 
when t > t 



it is easy to verify that 7(1, v, t) = g^ 1 . 

The second part of the theorem states that any reachable state can be reached in 
arbitrary short time interval if we are able to use the class of unrestricted controls. 
Assume that g — 7(1, u, t). For any given time t\ > define a control v as 

v(r) = £•«(£), rtR+ 

It is easy to show that j(l,v,ti) — g which implies that 72(1, t) C 72(1, ti) for 
arbitrary t\ > 0. □ 
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Summarizing the essential properties of the homogeneous control systems we 
have the following 

• The set reachable from 1 is a subgroup of G; 

• The set reachable from 1 is one and the same for three classes of controls: 
every g £ G that can be reached from 1 by means of unrestricted control 
can also be reached by means of only "bang-bang" control (possibly later); 

• If we use the unrestricted controls then every reachable state g £ G can be 
reached in arbitrary short time interval. 

For the affine case the most important facts are that if the set reachable from 1 
is dense in the Lie group G corresponding to the Lie algebra generated by the 
set {Xo, . . . , X m } then 7^.(1) = G and if the subgroup G is compact then always 
K(l) = G and 3 t > : K(l, t) = 11(1) (see 0). 

A control system of the typeEJis said to be controllable from g 6 G if IZ(g) = G. 
It is said to be controllable if it is controllable from every g G G. 

For right-invariant control system on Lie group the controllability from 1 is 
equivalent to the controllability because as we know TZ(g) = TZ(1) ■ g. As it follows 
from the reachability criteria a necessary condition for the controllability is that 
the group G is connected and g — g. If G is compact or the system is homogeneous 
then this condition is sufficient too. 

Though the class of affine (or homogeneous) Hamiltonians is very restricted, 
actually the same methods can be used for the reachability problem when the 
Hamiltonian is a continuous mapping H : R m — > g and as it will be evident 
from the further discussions, for a compact Lie group more important is not the 
Hamiltonian itself, but the class of controls and the image set of the Hamiltonian. 

7. Control System on Compact Lie Group with Continuous 

Hamiltonian 

This section is mainly based on the results of the work : 8, the main idea of which 
is that if we weaken the condition for the Hamiltonian to be affine, or homogeneous, 
and consider any continuous Hamiltonian, for compact Lie groups many essential 
results are still true. More formally, consider a right-invariant control system 

g = H(u)g, g€G,u£C, g(0) = 1 

where G is a compact Lie group, g is its Lie algebra, H : R m — ► g is a continuous 
map, C is the set of controls consisting of all piecewise-continuous maps from R + to 
K m , and H(u)g denotes the tangent vector at g 6 G obtained from H(u) S g by the 
right action of g. Let g be the Lie subalgebra of g generated by Image(iJ) C g and 
G be the corresponding connected Lie subgroup in G. The latter is the maximal 
integral submanifold for the right-invariant differential system defined by g. Since 
the control system is right-invariant and each trajectory of this system is tangent 
to g at the point 1 6 G, we have that they are contained in G. Therefore, the set 
reachable from 1 is a subset of G. After this, without restriction of generality it 
can be assumed that g = g and G — G. 

Using the same method as for the case of affine Hamiltonian (see Theorem EJ, it 
is easy to show that the reachable set 1Z(l) is a semigroup. 

To proceed further we need the following 
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Theorem 11 (see [H])- Let G be a connected Lie group and X = {Ai, . . . ,X rn } 
be a set of generators of its Lie algebra g. Then every g £ G is a finite product of 
elements of the form exp(iAi), t S M, G X . 

Proof. The set of all finite products of the form exp(tA^) is a path-connected sub- 
group G' of G. According to the Yamabe's theorem (see EJ G' is a Lie subgroup 
of G. Since G' contains the elements of the type exp(iAj) its Lie algebra g 1 con- 
tains the set X, therefore it contains g. Consequently, we have that g' = g. As 
the subgroups G and G' are connected with one and the same Lie algebra, they 
coincide. □ 

Since g is the Lie algebra generated by Image(iJ), we can take the set of gen- 
erators X from Image(TJ). For the reachable set 1Z(1) the elements exp(tXi) are 
constructed by using of only positive time t, while for the elements of the group G 
can be used any tel. 

Lemma 12 (see [S]). The set reachable from 1 is dense subset of G. 

Proof. Consider an element g 6 G such that g = J}exp(ijA Pi ) where all ti-s are 

i 

negative. We have that h = g^ 1 is an element of 72(1) (because all times become 
positive). Consider the sequence h n , n € N. Since the group G is compact this 
sequence contains a convergent subsequence {h nk \ k £ N}. We can assume that the 
sequence of subscripts {nk} is increasing. Take the sequence {h' h — /i n *+ 1-n * -1 }. 
Since n& is assumed to be increasing, we have that n k +i — nu — 1 > and therefore 
{h' k } C 72.(1). But obviously lim h' k = h^ 1 = g. We obtain that the elements 

k — >oo 

g = Y[ exp(tiX Pi ) with only negative ti-s can be approximated by the elements of 

i 

72(1) which implies that 72(1) is dense in G. □ 

The following lemmas help us to go further and show that the set reachable from 
1 coincides with G. For g £ G and V £ g let Ad(g)V be the element of g obtained 
by the adjoint action of g. 

Lemma 13 (see 0). Let X — {Ai,...,X m } be a set of generators of the Lie 
algebra g. If m = dim(g) then for any e > there exist two elements X p and X q in 
X and a time r, \t\ < e such that the element X m+ i = Ad(exp(TX p ))X q is linearly 
independent from {Ai, . . . , X m }. 

Proof. From m < dim(g) follows that at least for one pair of elements X p and X q , 

their commutator [A p , Xq] is linearly independent from the set A. If we assume 

m 

that Ad(exp(tX p ))X q = ^ <Zj(i)Aj, then the derivation of this equality by t gives 

i=l 

m 

[Xp,X q ] = o,i(0)Xi which is a contradiction. □ 

4=1 

Applying the same procedure to the extended set of generators {Ai, . . . , X m , A m+ i 
and etc., we can formulate the following 

Lemma 14. The set of generators {Ai, . . . , A m } can be extended to a basis of g 
- {Ai, . . . , X m , A m +i, . . . , X m + p } so that each element X m +i, i = 1, . . . ,p can be 
written as 

X m+i = Ad(exp(t q Xi g ) ■ ■ ■ exp(t 1 X h ))Xi 
where {Xi , . . . , Xi q } is a subset of {Ai, . . . , X m }. 
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Moreover, in the work of D. D'Alessandro (see [H]) there is calculated some upper 
bound for the number of factors in exjp(tqXi g ) ■ ■ •exp(tiX; 1 ). 

Lemma 15. The set reachable from 1 contains an open subset of G. 

Proof. Let {X±, . . . ,X n } be a basis of q as a vector space. Consider the mapping 
F : K" — >G 

n 

F(ti,...,t n ) = Y[exp(tiXi), (ti,...,t n ) eM" 

z=l 

It follows from the inverse function theorem that F is a local diffeomorphism on a 
neighborhood of G R™. This implies that for sufficiently small 5 > and some 
e G (0,6), the open neighborhood of (8, ... ,8) <E R m - 

{(ti,...,t m ) \8-e<U<8 + e, i = l,...,ra}cM™ 

is diffeomorphic to some open neighborhood V C G of g = F(5, . . . , 8) G G. Ac- 
cording to Lemma ITU there is a basis of — {Xi, ■ ■ ■ ,X n }, where the first m 
elements are generators from Image(TJ) and after the subscript m the elements are 
of the form 

X m+k = g k X Pk g~ 1 , p k G {l,...,m}, g k G G, k=l,...,n-m 
Regarding this, we obtain that the mapping F is of the form 

F(ti, t n ) = e tlXl ■ ■ ■ e tmXm gi e tm + lX ^ g^g 2 e tm+2X ^ g? 1 ■ ■ ■ 

■ ■ ■ g k e tm+kXp " gfc^gfe+i • ' ■ 9n- m e tnXp "- m g^n 

Consider the following terms of this product: 

.9i = 9i, 9i = 9i 1 92, ■ ■ ■ ,g' k = 9k-i9k, ■ ■ ■ ,9n-m = 9n-m 

They are not necessarily in but since 7^.(1) is dense in G we can select elements 

of 1Z(1) - hi, ... , h n -m arbitrarily close to g[, . . . , g' n ^ m 

hi ~ gi, ■ ■ ■ , h n — m Ri 9n—m 

and consider the mapping F : M™ — > G 

F(h, ...,t n ) = e tlXl ■ ■ ■ e t ™ x ™/ ll e* m + lX " h 2 ■ ■ ■ e*" Xp "-™ h n - m 

The latter is close to F and therefore maps some open neighborhood of (8, ... , 8) 
onto an open neighborhood of g' — F(8, . . . ,8). But for positive ti-s we have that 
F(ti, . . . , t n ) is an element of 71(1). Hence, the set 7^.(1) contains an open subset 
of G. □ 

Lemma 16. // the set reachable from 1 contains some open subset of G then it 
contains an open neighborhood of 1 . 

Proof. Let V be an open subset of G contained in 1Z(1). Consider the set V~ x = 
{9 1 I 9 e V}- K is clear that V^ 1 is also an open subset in G, though it is not 
necessarily contained in 7^(1). Since 7\L(1) is dense in G (see IT2|) . the set V^ 1 
contains at least one element h G TZ(1). Consider the set h ■ V. As h G V^ 1 clearly 
the set h ■ V contains 1 and is open in G. Because h e 7Z(1) and V C 7\L(1) and 
TZ(1) is a semigroup, we have that h-V& 7^.(1). Therefore, TZ(l) contains an open 
neighborhood of 1. □ 
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Hence, we have that 72.(1) is a semigroup in G and contains an open neighborhood 
of 1 G G. These implies that 72.(1) = G. 

To summarize, we can conclude that in the case of continuous Hamiltonian 
and unrestricted controls on a compact Lie group, the set reachable from 1 is the 
connected Lie subgroup corresponding to the Lie algebra generated by Image(H) 
( regardless of the structure of the Hamiltonian H : R m — > g ). 

Remark 5. Actually, during the proofs of the lemmas we use only "bang-bang" con- 
trols which is in accordance with the fact that the reachable sets for the unrestricted 
and "bang-bang" controls are one and the same. 

8. Minimal Set of Generators for su(n) 

As it was discussed in the previous section the controllability for right-invariant 
"bang-bang" (or piecewise continuous) systems on compact Lie groups is completely 
described by the differential systems defined by the set of generators contained in 
the image of the corresponding Hamiltonian. Hence, the question about a minimal 
set of generators for a given Lie algebra (and the corresponding Lie group too) is 
natural. For this question the following classical results of Kuranishi (see ^0]) and 
more resent results from ^2], ^2] and 52 are quite useful. In this section 
we give a review of some of them. First let us recall the general structure of a 
semi-simple Lie group and the corresponding Lie algebra. 

Let G be a compact Lie group and g be its Lie algebra. The rank of G is 
defined as the dimension of the maximal abelian subalgebra t) of g. This number 
coincides with the dimension of the maximal torus (i.e., the maximal connected 
abelian subgroup) H in G. 

All irreducible representations of H are 1-dimensional and are described by the 
elements of the group of characters H = Hom(H, C*), where C* denotes the group 
C \ {0}. Let g c denote the complexification of g: g c — g ® C. Any representation 
R of G in a complex vector space V defines a representation R of g c in V . The 
restriction of R on H gives the decomposition V = (J) V\ , where A-s are some 

A 

elements of (f) c )* and V x = {v e V \ R(exp(h))v = exp(A(/i)) -v, h G f) c }. 

A vector A G (f) c )* is called a weight of the representation R if dim(Vx) > 0. 
The dimension of V\ is called the multiplicity of the weight A. 

The adjoint action of G on g c is a representation of G on g c the restriction of 
which on H gives the decomposition 

where a G (f) c )* and 

g Q -{XG0 C | [Y,X] = a(Y)-X, VFGf) c } 

The nonzero weights of the adjoint representation are called the roots of the group 
G with respect to the maximal torus H . Let A denotes the set of roots. 

Let ( , ) be a G-invariant scalar product on g. For instance, it can be the Killing 
form (X,Y) = tr(adx ° ady). This scalar product can be linearly extended to 
g c . The algebra g is said to be semi- simple if the scalar product ( , ) on g is 
non-degenerated. In this case the group G is said to be semi-simple too. Further 
we assume that the group G (and hence, the Lie algebra g) is semi-simple. 
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3a, ep] = 



Notice that all roots are pure imaginary: A C it)* and for any a £ A we have 
that g a = Q- a - The latter implies that A = —A. Besides that dim(g Q ) = 1. The 
following theorem summarizes the structure of a semi-simple Lie algebra. 

Theorem 17. Let G be a connected compact Lie group with dimension n and rank 
I. Then G has 2m roots {ia k | k = 1, . . . , m}, where n = I + 2m. The group G is 
semi-simple if and only if it has I linearly independent roots. If G is semi-simple 
we can choose vectors e a £ l a , a £ A so that 

[e Q ,e_ Q ] = h a € f) 

and 

if a + (3 is not a root 

Ca0 e a+f3 if oi + (3 is a root 

The set {h a \ a £ A} contains a basis for f) and the vectors e a are called the root 
vectors. 

An element h € F) is said to be regular if {X 6 g | [h, X] = 0} = f). Otherwise 
the element h is said to be singular. It is clear that h £ f) is regular if and only if 
a(/i) ^ for each a € A. 

An element h £ E} is said to be strongly regular if ai(/i) ^ ct2{h) for any two 
distinct roots «i and «2- Prom the property A = —A follows that if h is strongly 
regular then it is regular too. 

Based on these data about the structure of a semi-simple Lie algebra we present 
the proof of the following classical theorem (see 

Theorem 18. Let L be a semi-simple Lie algebra. There exist two elements in L c 
which generate the entire Lie algebra L c . 

Proof. Let A be a maximal abelian subalgebra of L and L c = A c + l a be the 

ceGA 

corresponding root space decomposition of L. Let one element we are looking for 
be any strongly regular element h £ A and the other one be the vector e = ^ e a . 

Consider the vectors 

e = e 

ei = [h, e ] = a(h)e a 



e k+ i = [h, e k ] = a(h) k e c 



e2m-l = [K <Z2m-2[ = J2 a(h) 2m 1 e a 

where 2m is the number of elements in A. Since the numbers a(h), a £ A are 
different, the matrix 



/I ... 1 \ 

ai(h) ••• a 2m (h) 
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where {ai, . . . , ct2m} = A, is non-degenerate and therefore, the set of vectors 
{eo, . . • , e2m-i} is a basis for (J) Z Q . Thus, from the vectors eo, . . . , e2 TO -i wc 

ctGA 

can construct the root vectors e a , a <E A and then use the vectors [e Q , e_ a ], a G A 
to construct the basis for A c . □ 

Notice that the above theorem concerns not the Lie algebra L itself, but its 
complexification. Now we consider the ways for constructing the pair of generators 
for a Lie algebra itself. We concentrate on the Lie algebra su(n). First we review 
the structure of its complexification in the framework of root space decomposition. 

Let A be any element of su(n) . We can take A as a diagonal matrix, otherwise 

n 

we can diagonalize it. So, we have A — diag(X\, . . . , X n ), ^2 A? = 0, Xj S i ■ R. 

i=l 

Let H be the maximal abelian subalgebra of su(n) c containing A, and let H' be 
the zero eigenspace of ad^: 

W = {X esu{nf | [A,X] = 0} 

It is clear that 

H = H ■<=>• \ p ^ X q , when p ^ q 
Therefore, the element A is regular if and only if its eigenvalues Xj, j — 1, . . . ,n 
are distinct. We assume that the latter is true for A and thus H is the subgroup 
of all diagonal (with respect to the basis in which A is diagonal) traceless matrices. 
We have n linear functions Xj : H — > C, j = 1, . . . , n, where for any X £ H we 
have X — diag(Xi(X) 1 . . . , X n (X)). Let £ pq be the matrix with (p, g)-entry equal 
to 1 and other entries equal to 0. It is clear that for any X € H we have [X, £ pq ] = 
a pq (X) -£ pq , where a pq (X) = X p (X) — X q (X). Therefore, the matrices £ pq € su(?i) c 
can be considered as the root vectors and the functionals a pq = X p — X q as the roots. 
Each root vector has the multiplicity 1 if and only if a pq ^ ai s for (p,q) ^ (I, s). 
It is so when H contains such matrix B that X p (B) — X q {B) ^ Xi(B) — X S (B) for 
(p, q) ^ (I, s). Recall that such element B is said to be strongly regular. We assume 
that the matrix A itself is strongly regular. 

The roots ai2, 023, . . . , a n -in are called the fundamental roots since the others 
are obtained as sums of pairs of them. The basis for su(n) c corresponding to the 
roots is 

{E k = £ kk - £k+ik+i I k = 1, , . . . , n - 1} U {£ pq | p, q = 1, . . . , n, p^q} 

This basis is known as the Wayl basis for su(n) c . It is clear that the elements of 
this basis are not in su(n) but such basis can be constructed from the Wayl basis 
as 

{h k = i ■ E k } U {U pq = £ pq - £ qp } U {V pq = i ■ (£ pq + £ qp )} 
It is worth to notice the following commutating relations: for any diagonal A — 
i ■ h 6 su(n) we have 

[A, Upq] = a pq {A)V pq , [A, V pg ] = -a pq (A)U pq 

(4) 

[Upp+i, Vpp+i] = 2h p 

Our further discussion in this section is in accordance with |12j . Hence, we have a 
matrix A 6 su(n) which is a regular element (i.e., its eigenvalues have multiplicity 
1) and therefore, the associated Cartan subalgebra in su(n) is the set of all traceless 
diagonal matrices. Let B be another element of su(n). Consider the graph Ta,b 
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defined as follows: Ta,b has n nodes {|1), . . . , \n)} and the nodes \p) and \q) are 
joined if and only if the (p, g)-entry of the matrix B is nonzero. The graph is said 
to be connected if for all pairs of nodes there exists an oriented path connecting 
them. The following two theorems are very important for the controllability of 
affine systems on su(n) (see |12j 1 

Theorem 19. For A, B £ su(n) if A is diagonal, a necessary condition that the 
pair {A, B) generates the entire su(n) is that the graph Ta.b is connected. 

Theorem 20. Given a pair of matrices A, B £ su(n), assume that the graph Ta.b 
is connected. If A is strongly regular then the pair (A, B) generates the entire 
Lie algebra su(n). 

The first step of the proof (as in the case of the Theorem ll8|l is to consider the set 
of matrices A, [A,B],...,[A,..., [A, B]] and then use the commutating relations |3| 



n 2 ~ n— 1 

to construct the remaining n — 1 elements of the basis by the commutators between 
them. 

9. Quantum Gates: Recursive Construction of Generators 

Regardless of these, in quantum computation, it is important to construct the 
set of generators (gates) recursively: assuming that we have the complete control- 
lability of some subsystem of a large system and the task is to construct an optimal 
extension of the controllability to the entire system. Our further discussions makes 
more exact this question. 

Let G be a compact Lie group with Lie algebra g; go be a Lie subalgebra of g 
and Go be the corresponding connected Lie subgroup in G. Assume that we have a 
set of generators of go, thus the set reachable from 1 £ G by using these generators 
is Go- The problem is to find a minimal extension of the given set of generators 
of go to a set of generators of g. This problem is actual for control systems used 
for quantum computations, because, usually systems of universal quantum gates 
are constructed hierarchically, by adding new gates to existing lower-dimensional 
subsystems. 

Consider the adjoint action of the subgroup Go on the Lie algebra g. Since go 
is the Lie algebra of Go, it is invariant under this action and thus we have the 
following decomposition g = go © L\ © • • • © L p where each Li is invariant and 
irreducible for the adjoint action of Go. 

Remark 6. For instance, we can take the orthogonal complement of go under a 
scalar product invariant for the adjoint action. This space itself is invariant for the 
adjoint action of Go, and then take its decomposition on irreducible components. 

Lemma 21. Consider a set A = go U {Xi, . . . ,X p } where Xj £ Li \ {0} , i = 
l,...,p. The set A generates the entire Lie algebra g. 

Proof. Let g' be the Lie algebra generated by the set A. Since A is a subset of 
g', the space g' is invariant under the adjoint action of the subgroup Go- For 
each i = 1, . . . ,p the intersection g' n Li is nontrivial, because it contains at least 
Xi. Since the subspaces g' and Li are invariant for the adjoint action of Go, their 
intersection is also invariant. But by assumption the adjoint action of Go on each 
Li is irreducible, which implies that Li is a subspace of g' for each i = 1, . . . ,p and 
consequently g' = g. □ 
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After this it is still not clear is the set {Xi, . . . ,X P } a minimal complement of 
go to a set of generators of g or not. It is so when p = 1, but in general the answer 
is negative. To illustrate this, consider the case when dim(Li) = 1, i = 1, .. . ,p. 

Let us select nonzero vectors e.; G L^, i = l,...,p. The set {ei,...,e p } is a 

p 

basis for L = @ Li, and we have that [x,ei] = Xi(x) ■ e^, x G go, where Xi is 

i=l 

the character of the representation of Go on Li via the adjoint action. Consider 

p 

the vector e = ^2 ei. If we are able to find such elements Xi, . . . ,x p in go that 

i=l 

the matrix I a\ — Xi(xj) ) is non-degenerate then the vectors [xi, e], . . . , [x p , e] 

form a basis for L and therefore we need just one element e which together with go 
gives a set of generators of g, and of course it is a minimal complement for any set 
of generators of go to a set of generators of g. 

The method described in the Lemma 1211 is in agreement with the standard 
method of constructing the universal set of gates for quantum computation. Now 
we describe a method of recursive construction of generators follow the ideas dis- 
cussed in Let X and Y be finite-dimensional Hilbert spaces. Assuming that 
we have universal sets of generators for the unitary groups U(X) and U(Y), the 
problem is to extend them to the universal set of generators for the unitary group 
U(X <g) Y). First let us recall some useful facts from the representation theory. 

Let V be a finite-dimensional Hilbert space and tt : G — ► U(V) be a unitary 
representation of a Lie group G. For g G G and v G V we shall use the notation gv 
for the action w(g)(v). Let us denote by C(tt) the set of all such linear operators 
A : V — ► V that Aon(g) = ir(g)oA, V g G G. It is clear that C(tt) is a subalgebra 
of the algebra of all endomorphisms of V. 

Lemma 22. The representation n is irreducible if and only if the algebra C(tt) 
does not contain a projection operator different from and 1. 

Proof. If the representation tt is reducible then V = Vi © V2 where V% and V% 
are proper non-trivial invariant subspaces. It is clear that the projector on the 
subspace V\ (or V2) commutes with every n(g), g G G. Conversely, if P : V — ► V 
is a non-trivial projector on a proper subspace which commutes with all 7r(g), then 
we have 

V x G Im(P), V g G G : P{gx) = gP(x) = gx =>■ gx G Im(P) 

Hence, the subspace Im(P) is invariant for every n(g) and therefore the represen- 
tation 7r is not irreducible. □ 

Lemma 23. The representation n is irreducible if the algebra C(tt) does not contain 
a self-adjoint operator different from and 1. 

Proof. If 7r is not irreducible then according to the previous lemma the algebra 
C(tt) contains a projection different from and 1, which is a self-adjoint opera- 
tor. Conversely, if C(vr) contains a self-adjoint operator A different from and 1, 
consider the spectral decomposition of A: A = XiPi, where Pj-s are orthogonal 
projectors. From Ag = gA follows that if x G Im(Pi) then we have 

A(gx) = gA(x) = Xi • gx => gx G Im(p) 

This implies that each subspace Im(Pi) is invariant for the representation tt. □ 
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The last lemma implies the following 

Theorem 24. The representation it is irreducible if and only if the algebra C(n) 
consists of only scalar operators (i.e., operators of the form A • 1, A 6 C). 

Proof. Let A S C(ir) be a non-scalar operator. Then at least one of the following 
two self-adjoint operators B\ = A + A* and £? 2 = i(A — A*) is also non-scalar. 
According to the previous lemma this implies that the representation 7r is not 
irreducible which contradicts to the assumption of the theorem. □ 

For two representations %i : G\ — ► U(V\) and tt 2 : G 2 — ► ^(^2)1 where G\ 
and G 2 are Lie groups and V% and V2 are Hilbert spaces, let us denote by 7Ti ® 7r 2 
the representation of G\ x G2 in V\ ® V2 defined as 

(tti ®TT2){gi,g2){x,y) = g x x®g 2 y, gi € Gi, g 2 £ G 2 , x € Vy, y 6 V 2 

Theorem 25. J/ the spaces V\ and V2 are finite- dimensional then C{"K\ ® 1x2) — 
C(7ri)®C(7r 2 ). 

p 

Proof. Any operator T : V\ <8> V2 — > V\ <8> V2 can be represented as T = A p Cg) B p 

where A p is an endomorphism of V\ and B p is an endomorphism of V 2 . The RH 
side of this equality can be reduced to the form when the operators B p are linearly 
independent. If T £ C(7Ti ® 7r 2 ) then we have the following 

[T,gi ® 52] =0, V gi e Gi and V g 2 e G 2 

This implies that 

[T,0®1]=O, V g e Gi 

Expanding T we obtain 

p 

®B P = 0, V.geGi 

Since the operators £? p are linearly independent, we have that each [A p ,g] is 0. 
Therefore, each A p S G(7Ti), which implies that T € End(V\) ® G(7T 2 ). In same 
manner we obtain that T G G(7Ti) ® End(V2). Thus, we have that T e C(tti) ® 
End{V 2 )^End{Vi) ® G(tt 2 ) and therefore T e G(tti) f| G(tt 2 ). □ 

Corollary 1. if the representations %\ and n 2 are irreducible then the representa- 
tion 7Ti ® 7T2 is ako irreducible. 

Now consider the unitary group U (Vi^Va) where Vi and V2 are finite-dimensional 
Hilbert spaces and dim(Vi) = to, dim^) = n. First of all notice that its Lie algebra 
is 

u(Vi ® V2) = u(Vi) ® iu(V 2 ) 

because (see US]), obviously the Lie algebra u(Vi) £S> iu(Va) with the commutator 
defined as 

[01 <S> ibi,a 2 <£> ib 2 ] = a 2 ai ® 6 2 6i — aia 2 (8) 61 6 2 
is a subalgebra of u(Vi ® V2). Also we have that 

dim(u(Vi <g> y 2 )) = (tow) 2 = dim(yi) 2 dim(V 2 ) 2 = dim(u(Vi) ® ni(V 2 )) 
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Consider the subgroup of U(Vi © V2) generated by the elements of the form g\ © 
.92, 51 E U(Vi), 32 G J7(F 2 ). Let us denote this subgroup by U(V\) © f7(F 2 ). The 
Lie algebra of this subgroup is the subalgebra of elements of the form 

u± © 1 + 1 © m 2 , ui £ u(Vi) u 2 € u(V 2 ) 

Such element is known as the Kronecker sum of u\ and w 2 . 

Remark 7. The Kronecker sum oju\ andu 2 can be written as ui©z(— il)— il®iu 2 
in accordance with the equality u(Vi © V2) = u(Vi) © 111(^2). 

Consider the decomposition 

u(Vi © V2) = (iK • 1 © fiu(Vi)) © (R • 1 © i ■ su(V 2 )) = 

= (iR® 1) © (su(Fi) © 1) © (1 ©su(F 2 )) © (su(Fi) © i ■ su(V 2 )) = 

= (u(Vi) © 1 + 1 © u(V 2 )) © (su(Vi) © i ■ su(V 2 )) 

The first summand is obviously the Lie algebra of U(Vi) © U(y 2 ) = G. Let us 
denote this Lie algebra by L. Consider the adjoint action of the group G on the Lie 
algebra u(Vi © V 2 ). It is clear that L is invariant under this action. Since the adjoint 
action of U{V\) on Su(Vi) is irreducible and the same is true for U (V2) and i •su(V 2 ), 
we have according to the corollary ^ that the action of G on su(Vi) © i ■ su(V 2 ) is 
irreducible too. This together with Lemma |^ implies that if we have a set of 
generators of L we need only one new element from su(Vi) ©i -suV^ to generate the 
entire Lie algebra u(Vi © V2) (actually, any one element outside of L is enough). 
Hence an optimal (in the sense of quantity) extension of the set of generators of 
the Lie algebra of U(Vi) © £7(14) to the set of generators of the Lie algebra of 
U(Vi © V 2 ) is obtained by adding one and only one element outside of the Lie 
algebra of U(Vi) ® U{V 2 ). 

The set of generators of Lie algebra is used to construct a control system the cor- 
responding set reachable from 1 of which is the entire Lie group. But the dynamics 
of a control system is a continuous process while the conventional computational 
process (quantum too) is discrete. In other words, the latter uses not 1-parameter 
flows of the form exp(tXi), t > as gates to generate group elements, but discrete 
"flows" of the form gf, n = 1, 2, . . ., where gi-s are so called gates. The generators 
of Lie algebra can be used to construct such gates. 

Let X = {Xi, . . . , X m } be a set of generators of the Lie algebra L of a compact 
connected Lie group G. For any Xi £ X consider the corresponding 1-parameter 
flow exp(tXi), tel. If exp(tXi) is periodic and therefore exp(TXi) = 1 for some 
T G M + then take the element gi = exp(aiTXi), where cti G K + is an irrational 
number. The set S(gi) = {g™ \ n £ N + } is dense in the orbit exp(tXi), t G M. 
If for some generator Xj £ X the corresponding orbit exp(tXj), t £ IR + is not 
periodic then consider the maximal torus C G where it is contained, replace 
the element Xj in the set of generators X by periodic generators of the Lie algebra 
of the torus T fe . The obtained set is again a set of generators of L because Xj is 
a linear combination of the generators of the torus. In such way we can achieve 
that all the generators of L are periodic, and then construct elements of the group 
G: gi,...,g m , such that for each of them the group S(gk) = {g% I n £ N + } , k = 
1, . . . , m is dense in the 1-parameter subgroup exp(iXfe), tel. Since the subgroups 
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exp(tAfc), k — 1, . . . , m generate the group G, the subgroups S(gk), k = 1, . . . , m 
generate a dense subgroup in G. 

Lemma 26. If G is a Lie group with representation in some Hilbert space TL and 

V C 7i is a subspace invariant for some dense subset S C G then V is invariant for 
the entire group G. 

Proof. Consider the continuous map F : G — ► Hom(V, V ) defined as F(g) = 
P o g\ v where P : H — > V 1 - is the operator of orthogonal projection on V . Since 

V is invariant for S, we have that F(S) — 0. Since S is dense in G and F is 
continuous, we have that F(G) = which means that V is invariant for the entire 
group G. □ 

Let L\ be a Lie subalgebra of L and G\ C G be the corresponding connected Lie 
group. Let L — L\ © A be such decomposition of L that the subspace A is invariant 
under the adjoint action of G\ and the action of G\ on A is irreducible. As it follows 
from Lemma l^Tl if we take only one element X £ A, then the set Li n{A} generates 
the entire Lie algebra L. If X is such element that the 1-parameter subgroup 
{exp(iA) \ t eR} is periodic, then for some a G R, the element g = exp(aA) 
generates a dense subgroup in {exp(iA) | t e R}. In this situation the set G\ n {g} 
generates a dense subgroup in G. Moreover we the following 

Theorem 27. Under the above conditions the set G\ (~l {g} generates the entire 
group G. 

Proof. First let us show that if go = g m is a normalizer of the subgroup G\ then 
m = 0. Assume that m > and goGig^ 1 — G\. It is clear that in this case the 
subgroup {#q n G Z} is also a normalizer of Gi. But this subgroup is dense too 
in {exp(tX) \ t € M} = S. Consider the adjoint action of S on the Lie algebra L. 
Since {g^ \ n £ Z} is a normalizer of Gi, we have that the group S contains a dense 
subgroup for which L\ is invariant. This implies according to the previous lemma 
that L\ is invariant for the adjoint action of the group S. Therefore we have that 
[A, L\] C L\. Since A is invariant under the adjoint action of G\ and X £ A, we 
have that [X, L%] C A. The result is that [A, L{\ = 0, which contradicts to the 
assumption that the action of G\ on A is irreducible. Hence, we obtain that there 
is no nontrivial normalizer for G\ in {g n \ n £ Z}. 

Now take any element g\ = g m , m/0 and consider the group G2 = giGig^ . 
Since g\ cannot be a normalizer of Gi, we have that Gi ^ G\. Therefore, the 
group G generated by Gi H G2 contains Gi as a proper subgroup. The group G is 
connected because G\ is such. Therefore, its Lie algebra L contains L\ as a proper 
subalgebra. From this follows that L n A ^ 0. But as the adjoint action of G\ on 
A is irreducible and L is invariant under this action, we have that L fl A — A. The 
latter implies that L = L. Since G and its subgroup G are connected and they have 
a common Lie algebra, they must be equal: G = G. □ 

Regarding the above theorem in the context of inductive construction of genera- 
tors of the unitary group of tensor product, we have the following situation: consider 
the unitary group U(Vi ® •• • <8 V m ), where Vi,...,V m are some finite-dimensional 
Hilbert spaces; if we are able to generate the subgroups U(Vi), . . . , U(V m ), then we 
can construct a set of generators of the entire group step-by-step by adding only 
one element at each step. For instance, consider the group U(V\ <g> V2 <8> V3 ® V4 ® V5. 
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If we have the groups U(Vi), U{V 2 ), U{V 3 ), U{V 4 ) and U(V 5 ), then to construct 
U(Vi ® V 2 ) and U(V 3 ® 14) we need only 2 more elements. Then with U(Vi <8 F2) 
and t/(V3 £g> V4) we can construct U(Vi <8> V2 ® V3 <& V4) by adding 1 more element. 
And finally, with U(Vi C5> V 2 <8> V3 (E> V4) and 17 (V5) we can construct the entire group 
U(Vi <8> V2 ® V3 ® V4 ® V5) by adding again 1 more element. Therefore, besides 
the subgroups U(Vi), U(V 2 ), U{y 3 ), U(V/t) and U{V$) we need 4 more elements 
to construct the group U(V\ <g> ® V3 <g) V4 ® V5). 

The above construction is, of course, in accordance with |15|. and the key point 
here is that an element outside of the subgroup is not a normalizer of this subgroup. 

10. Control Systems on Homogeneous Spaces and Grassmann 

Manifold 

Let M be a smooth manifold and G be a Lie group acting on M. Let g be the 
Lie algebra of G. Any element X £ g generates a vector field tt(X) on M defined 
as follows: for m G M let 

7T(X) m = R' m (l)X 

where R m : G — > M is the map R m (g) — gm, g 6 G. We shall use the notation 
Xm for the value of the vector field n(X) at a point m G M. The vector field 7r(X) 
is not invariant under the action of G and its conversion rule under this action is 
the following: for g € G we have 

g(Xm) = gXg-'gm, VmeM =>■ s(tt(X)) = ^Ag- 1 ) 

We consider control systems on M of the form 

(5) m = 7r(A M ) m , m e M 

Where X u G g and u is a control parameter. Any solution of the system [S] is 
contained in one orbit of the action of G on M . This fact easily follows from the 
definition of the vector field tt(X): for any point m G M the vector 7r(A) m is 
tangent to the orbit of the point to. Though the vector field tt(X) is not invariant 
under the action of G but we have the following 

Lemma 28. The integral manifolds of the differential system on M generated by 
the mapping §9Xh tt(A) are exactly the orbits of the action of G. 

Proof. We apply the theorem of Nagano-Sussmann (see OJl to this case. For any 
X G the action of the 1-parameter flow exp(tX), t G K on the manifold M 
generates the vector field which is exactly the same as the vector field 7r(X). For 
any Y G g, by the conversion rule of the vector field 7r(Y) under the action of G we 
obtain: 

exp(<A)7r(y) = 7r(exp(iA)r cxp(-iA)) 

which implies that ir([X, Y]) = [tt(X), 7r(Y)]. Hence, the differential system m 1— > 
""(fl)m> to G M is closed under the bracket operation. The 1-parameter flows of the 
vector fields 7r(A), i£g (as it was mentioned just before) are m <— > exp(iA)TO, m G 
M, and again according to the conversion rule we have that for any g G G: 

gir{g) = ir{ggg~ l ) = tt(L) 

Therefore, the differential system 7r(fl) is invariant under the action of G, which 
implies according to the theorem of Nagano-Sussmann that the orbits of the action 
of G are the integral manifolds of the differential system ir(g). □ 
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From these discussions easily follows that a control system on homogeneous 
space has one significant property: the complete controllability of such system is 
local (moreover: it is pointwise): if G acts on M and the set of vector fields n(V) 
generated by some V C Q is completely nonholonomic in a point too e M, then 
tt(V) is completely nonholonomic in any point m e M. 

Now we consider the problem of controllability of finite-dimensional quantum 
systems in the framework of control systems on homogeneous spaces. The classical 
notions of transformation and state controllability can be carried to such systems. 
In this case the corresponding Lie group is the unitary group U(n) and the control 
system X u is a parameter-dependent element of the Lie algebra u(n) . 

The notion of state in quantum (as well as in classical) case has various meanings. 
The most general of them is density matrix, i.e., a positive self-adjoint matrix with 
trace equal to 1 . The equation for the control system on the set of density matrices 
is of the form 

P = [X u , p] 

where p is a density matrix. Since any solution of such system belongs to one 
orbit of the adjoint action of U(n): p(t) <E {gpog^ 1 I S 6 U(n)}, let us recall the 
structure of an orbit of density matrix. For any density matrix p consider its 
spectral decomposition 

k k 
p = ^A,P 4 , A 4 >0, J! A, = 1, P*=P h P?=Pi, i = l,...,k 

i=l i=l 

The spectrum of the density matrix (operator) p is the set of pairs 

Spec(p) = {(A i; TOj = rank(Pi)) | i = 1, . . . , k} 

Two density matrices p\ and p 2 are in one and the same orbit of adjoint action 
if and only if Spec(pi) = Spec(p 2 ). Clearly, it is more reasonable to consider the 
controllability not for the entire set of density matrices but for each particular orbit 
of the adjoint action of U(n). 

Definition 7. A system X u e u(n) is said to be state controllable for a given orbit 
S of the adjoint action of U(n) on the set of density matrices, if for any pair of 
density matrices po, pi € S there exists an admissible control u and a time T > 0, 
such that the solution of p = [X u , p], p(0) — po, satisfies p(T) = p\. 

We shall formulate some criteria of state controllability for the special case when 
the orbit of the adjoint action is the Grassmann manifold GVfc(C"). This is equiv- 
alent to the case when Ai = 1/k and rrii = 1 i = 1, . . . , k. For any linear operator 
A : C" — ► C™ let Ai 2 be the upper-right component of the decomposition 

\A 21 A 22 J 

corresponding to the decomposition C" = C k © C n ~ k . For a given subset & C u(n) 
let Lie(&) be the Lie subalgebra of u(n) generated by 6. 

Theorem 29. Let H : W n — ► u(n) be a continuous map and C(R + , R m ) be the 
space of piecewise continuous functions from M. + to R m . For a given Grassmann 
manifold 

Gr k (C n ) = {P:C > C | P* = P, P 2 = P, trace(P) = k} 
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the control system 

P = [H(u),P], P G Gr k (C n ), u g C(R+, R m ) 
is controllable if and only if for the subspace C fe — ► C n the set of linear maps 

X12 : C fe — ► C"" fe , X G Lie{Im(H)) 
gives the entire space Hom(C , C n ~ k ). 

The proof follows from the already mentioned fact that the necessary and suf- 
ficient condition for the complete controllability of a control system on a homoge- 
neous space is that the corresponding differential system be completely nonholo- 
nomic in at least one point, and the fact that the tangent space of the Grassmann 
manifold Gr&(C n ) at a point X is canonically isomorphic to Hom(X, X- 1 ). 
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